Tunneling through triple quantum dots with mirror symmetry 
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Indirect exchange interaction between itinerant electrons and nano-structures with non-trivial 
geometrical configurations manifests a plethora of unexpected results. These configurations can be 
realized either in quantum dots with several potential valleys or in real complex molecules with 
strong correlations. Here we demonstrate that the Kondo effect may be suppressed under certain 
conditions in triple quantum dots with mirror symmetry at odd electron occupation. First, we show 
that the indirect exchange has ferromagnetic sign in the ground state of triple quantum dot in a 
two-terminal cross geometry for electron occupation TV = 3. Second, we show that for electron 
occupation A'^ = 1 in three-terminal fork geometry the zero-bias anomaly in the tunnel conductance 
is absent (despite the presence of Kondo screening) due to special symmetry of the dot wave function. 
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I. INTRODUCTION 

Many-particle effects in quantum tunneling through 
quantum dots are extensively discussed in the current 
literature (see, e.g., recent reviewai^). Single-electron 
tunneling under the conditions of strong Coulomb block- 
ade is accompanied by co-tunneling processes with spin- 
reversal, which involve dynamical screening effects sim- 
ilar to the celebrated Kondo-scattering in magnetically 
doped metals. The pertinent effect was predicted and 
observed in single-well quantum dots with odd electron 
occupation, where the electrons confined in the well are 
represented by the non-compensated spin 1/2 of an elec- 
tron on the highest occupied discrete level. This generic 
pattern may be enriched in many ways, in particular, by 
studying the tunneling through complex quantum dots 
containing two or three valleys. 

In this paper we are interested in specific physical 
properties of electron tunneling through complex quan- 
tum dots containing several potential valleys with es- 
sentially different capacitances. Originally, the idea of 
coupling several nanoobjects having strong and weak 
Coulomb interactions is formulated in the context of elec- 
tron and spin structure of complex molecules (e.g., lan- 
tanocenes, containing strongly correlated f-electrons hy- 
bridized with weakly correlated molecular orbitals occu- 
pied by p-electrons). It is noticed^ that the energy differ- 
ence between the singlet (S) ground state and triplet (T) 
excited state of a molecule with even number of electrons 
N = Nf-\-Np is anomalously small, so that the triplet ex- 
citation affects the magnetic response of the system. The 
simplest artificial analog of this system is an asymmetric 
double quantum dot (composed of big and small dots) 



with even occupation (e.g., = 2) and charging ener- 
gies Qs ^ Qb for small and big dots. When this double 
quantum dot is coupled with the metallic leads via the 
big dot (the T-shaped geometry) then lead-dot electron 
tunneling may induce an S' — > T crossover^: the energy 
Ex of the triplet state becomes lower than the energy 
Eg of the singlet state, and the Kondo-regime may show 
up. In case of odd occupation = 1, the Kondo- Fano 
regime is relevant, where the Kondo effect induced by an 
electron localized in the small side dot affects the elec- 
tron tunneling through the weakly correlated big doA 
If the leads are connected via the small dot, the big dot 
plays the role of additional reservoir for Kondo screening 
in case of odd A^ and the two-channel Kondo effect may 
be realized under certain conditions&i 

The more complicated model of asymmetric triple 
quantum dot (TQD) with a small dot sandwiched be- 
tween two big dots was considered in Ref.^ for odd and 
even occupations A^ = 3, 4. It was shown that the Kondo 
regime is accessible in both cases. Moreover, the TQD 
with even occupation demonstrates SO{n) symmetry. 

Recently, a mirror symmetric TQD (see figure^ in an 
"open" regime (near the Coulomb blockade peak) were 
studied both experimentallj"^ and theoretically In this 
case the big "mesoscopic" dot is connected with two small 
dots and with metallic reservoirs. As a result, an indi- 
rect RKKY exchange interaction between the spins in 
the couple of small dots occurs, and its sign may be con- 
trolled by changing the parameters of the large central 
dot by applying an external magnetic field. 

Here we focus on a mirror symmetric TQD in the 
"closed" regime, that is, the valley between Coulomb 
blockade peaks where the total number of electrons in the 
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FIG. 1: Triple quantum dots in "cross" (a) and "fork" (b) 
geometry. 



TQD is fixed. Our main results are: (i) in the cross ge- 
ometry (Fig. ^) with electron occupation = 3, such a 
dot possesses an unusual property: the indirect exchange 
tunneling constant between the big dot and the leads is 
ferromagnetic; the Kondo tunneling is therefore absent 
although the TQD behaves as a local moment; (ii) in a 
"fork" configuration (Fig. ^) with electron occupation 
A'^ = 1, TQD exhibits two diff'erent tunneling regimes: 
depending on the gate voltages, the Kondo regime may 
be observed as a zero bias anomaly or as a finite bias 
anomaly in the tunnel conductance. 



II. TRIPLE QUANTUM DOT IN A CROSS 
GEOMETRY 

The TQD in the cross geometry (figure^) is composed 
of left I, center c, and right r dots, with corresponding 
levels and charging energies e.j,Qj,j — l,c,r. It is mod- 
eled by the Anderson Hamiltonian 



H = H, 
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containing the terms describing the dot, two leads and 
the dot-lead tunneling, respectively. The first term is 
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The parameters Qj are chosen in such a way that for elec- 
tron occupation = 3, and in the absence of inter-dot 
tunneling {W = 0), each dot is occupied by one elec- 
tron. In the mirror-symmetric case such configuration is 
realizable provided 

ei= er = es, A = ec - Es > 0, 

Ql ^Qr = Qs'> Qc ■ (3) 

At finite W charge transfer from the central dot to the 
side dots is possible, but the double occupation of the 
side valleys is still suppressed by strong Coulomb block- 
ade Qs ^ W. In the charge sector N = 3 and for the 



mirror symmetric configuration the lowest energy states 
are two spin doublets, (even and odd relative to the I <-!■ r 
permutation), a spin quartet state, and a doubly degen- 
erate charge transfer exciton (with two electrons in the 
central dot). The corresponding energies are 
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Here the inequality W//S. <C 1 is assumed to be valid. 
The eigenfunctions of the doublet and quartet states 
(which predetermine the structure of the effective ex- 
change Hamiltonian, see below) are: 
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5ij)\l^/2, and rotation 



Here b\^ = KV]^ - ^^^^(l 

angles 6'„ = arcsin (V3VF/A), 6g = arcsin(W/A). 

The two other terms in Eq. are the band Hamilto- 
nian describing the electrons in the leads 
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and the tunneling Hamiltonian 
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Following the standard procedure, one derives an indi- 
rect exchange interaction between lead and dot electrons 
by means of the Schrieffer- Wolff (SW) transformation. 
For a composite quantum dot, this is accomplished in 
terms of spin eigenstates of defined in Eqs. Q,©. 

Thus in the SW regime the main contribution to the 
lead-dot tunneling is given by the components ^ cos Og^u- 
The quartet state has the usual structure prescribed by 
a standard Young tableau for three electrons, {{ijk) in- 
dicates the cyclic permutation of three sites Icr) . 

After performing the SW-like canonical transforma- 
tion, one gets the effective exchange Hamiltonian 
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where Su is the spin 1/2 vector operator with compo- 
nents = \u ^){u i\,S^ = {\u T)(u T | - |zi i){u i |)/2, 
and s is the spin operator of lead electrons defined as 
s = X)/c '^efco-'^'^e/co-' ; T is the sct of Pauli matrices, Cea is 
the even combination of lead electron annihilation opera- 
tors (the odd one is excluded from the tunneling Hamilto- 
nian by standard rotation^^). Remarkably, the exchange 
constant J„ has ferromagnetic sign. 
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(led is the position of the central dot level relative to the 
Fermi energy of the leads, source and drain contacts are 
assumed to be equivalent, = = V). The reason 
for an unconventional sign of the exchange interaction is 
that only one of three electrons in the TQD is involved in 
exchange interaction with the leads, and the overlap of 
the two other electrons wave function entering the state 
I-Ducr) gives the factor —1. Thus in this geometry we 
encounter a unique situation, where the Kondo screening 
is ineffective for a quantum dot with odd occupation. 

Yet, a crossover from ferromagnetic to anti- 
ferromagnetic scenario is feasible. Indeed, the level spac- 
ing in the spin multiplet I0J is governed by the parame- 
ter W/ A. If this spacing is small enough, the renormal- 
ization of these levels due to lead-dot co-tunneling be- 
comes relevant in the RG flow equations along with the 
exchange screening in the framework of Haldane renor- 
malization procedurei^. Accordingly, the energy levels 
(0J) are renormalized as a result of integrating out the 
band edges and shrinking the band width from its bare 
value Dq to a smaller value D comparable with |ec|). The 
corresponding RG invariant is. 
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with tunneling rates T o(^u,g) ~ '"'/Oo cos^ and 
Tq « TrpoV^^- Due to the hierarchy Tq > > Tdu, 
a level crossing is feasible (see Fig. Physically, it 

implies a crossover from the non-Kondo (ferromagnetic 
exchange) regime to the under-screened Kondo regime 
with a pronounced maximum of the conductance around 
the degeneracy point. The parameters W and A, which 
determine the initial conditions for the flow equations 
(|1U|I can be controlled by gate voltages. Varying these 
initial conditions, one may tune the region of crossover 
to the SW regime at the point 
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For D < D the properties of the system are determined 
by the SW Hamiltonian. The effective Hamiltonian ^ 
may be written for D = (marked on the horizontal 
axis of Fig. If the condition Hll|l is fulfilled in the 
vicinity of the crossing point {D = Dcr), the exchange 
Hamiltonian (instead of JSJ) acquires the form 




FIG. 2: Haldane flow diagram for the levels -EaoI ljlO|l . d 
ttD/Tq. Energy is measured in meV units. 



expressed in terms of operators for localized spin 1/2, 
S„.g, spin 3/2, Sq, and the vector R which induces tran- 
sitions between the quartet \Q) and the doublet \Du)- 
There are no transitions between \Q) and \Dg) since these 
states have different I — r parity. 

At the point Dcr the degeneracy of spin state of 
the TQD is maximal, corresponding to the symmetry 
SU{2) X SU{2) X SU{2). Both to the right and to the left 
of this point some of the states in the spin multiplets are 
quenched at T ^ and depends on the energy gaps 
AQg = EQ{D)-EDg{D) and Ag^ = EogiD) - EDuiD). 
To find the function T/f (Aq^, Agu), one should solve the 
scaling equations for the coupling constants in the Hamil- 
tonian Hsw, 
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E[sw = Ju^u ■ s 



s + JqSq-s + JrR-s, (12) 



where ja = paJa (a = 9, u, Q, R), and po is the density of 
states which is assumed to be constant. The procedure is 
self-consistent because Tk itself predetermines the char- 
acteristic energy interval for states in the spin Hamil- 
tonian involved in its formation. Varying D in Fig. |21 
from the ferromagnetic non-Kondo regime {D ^ D^) to 
the crossing point 5^, one reaches the point > 
which arises due to influence of the excited states Eq{D) 
and Eug{D). Just then, Tk sharply increases, reaching 
its maximum value in the point of maximum degeneracy 
Dcr- Moving further to the left, the level Edu freezes out. 
This means that the vector R in the Hamiltonian H12|) 
does not contribute anymore to Kondo co-tunneling, and 
the Kondo effect is determined by the pair of states Eq 
and E^g, with the dynamical symmetry of TQD being 
SU{2) X SU{2). Further decrease of D eventually results 
in the quenching of Eog- The system then exhibits an 
under-screened Kondo regime of a localized spin 3/2 mo- 
ment. Fig. O illustrates these crossover effects on Tk- 



4 




FIG. 3: Evolution of Kondo temperature as determined by 
scaling equations 1)13^ . 



The evolution of Tk is reflected in the behavior of 
tunnel transparency and conductance as a function of 
energy. Far enough to the left and to the right of the 
crossing point, this behavior is stepwise. When the dot 
is in the doublet ground state {D ~ D^), the spin mul- 
tiplet as a whole contributes to the transparency at high 



energies uj ^ Aqg + Ag^- With decreasing u; the quar- 
tet Eq and the even doublet Ejjg energies freeze out in 
this order (they are not renormalized anymore). Even- 
tually, Kondo tunneling is quenched at low energies, so 
that the zero bias anomaly (ZBA) has a shape of a dip. 
On the other hand, in the crossover regime, the ZBA 
follows the conventional Kondo peak. Finally, the struc- 
ture of the peak at the regime D ~ Dq is even more 
complicated. Within the framework of our approach we 
may describe the evolution of transparency for T > 
where it is approximately described by the simple relation 
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FIG. 4: Tunnel transparency of TQD in a doublet (left panel) 
and quartet (right panel) ground states. T^eafc is the maxi- 
mum value of T{ui). 

This type of non-universal behavior of Tr- is known 
in the theory of strongly correlated quantum dots with 
even occupation, where the singlet-triplet level crossing 
usually occurs'*'^'^'^. The novelty of the present scenario 
is that it is manifested in quantum dot with odd occu- 
pation where the absence of Kondo effect occurs due to 
ferromagnetic exchange coupling with the localized spin 
doublet. The non- universality of Tx occurs as this ferro- 
magnetic exchange competes with two anti-ferromagnetic 



exchange interactions (with doublet and quartet local- 
ized moments), so that, in some sense, one deals with a 
"three-stage" Kondo effect. Thus, we have completed our 
discussion pertaining to the cross-shaped TQD of figure 



III. TRIPLE QUANTUM DOT IN A FORK 
GEOMETRY 



It is then natural to expect peculiar features of Kondo 
tunneling also in a mirror symmetric TQD in the fork 
geometry (Fig. ^) with Qi^r S> Qc- The dots and 
leads are labeled 1,2,3 and each dot is attached to its 
own lead. In case of = 3, the exchange coupling J3 
between the central dot and its adjacent lead is ferro- 
magnetic in accordance with Eq. whereas those for 
the two other (small) dots (Ji, J2) are antiferromagnetic. 
Besides, there are also non-diagonal exchange couplings 
Jij = Jji. All these are coupled within a system of RG 
flow equations. A question now arises, whether it is pos- 
sible to find a regime where the Kondo resonance arises 
only in dots 1,2, whereas dot 3 remains Kondo inactive? 

To answer this question, we calculate Tk and G within 
the same scheme as in the preceding section for a system 
described by the Hamiltonian with the tunneling term 
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instead of JT)). The effective exchange Hamiltonian of 
this system is obtained in the same way as (I12II . The 
mirror / — r symmetry entails Vi = V2 7^ V3. If the 
ground-state consists of the odd-parity doublet Ed^, , the 
corresponding SW Hamiltonian has the form 



Hsw = JjS ■ Si + J12S ■ (Si2 + S21) 

i=l 

+ J13S • (Si3 -I- S31) + J23S • (S23 + S32). 
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Here the exchange constant J3 < is the same as J„ 
whereas Ji_2 > 0. For cos6'„ « 1, these constants are. 
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The system of RG flow equations for the Hamiltonian 
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(|15|l now reads, 



dji 
dlnd 

din d 

dji2 
dlnd 
dji3 
dlnd 



-jiaUi +33)- 



(17) 



A somewhat unexpected results following the analysis 
of these equations: The coupling constants related to 
channel 3 become relevant in spite of their negative initial 
values The Kondo temperature is 



Tk = ^exp — 



(18) 
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where j+ = ji + ji2 + js, j- = ji + ji2 - h, so that 
the Kondo resonance arises in all non-diagonal channels 
and the TQD looses its "exotic" properties. The Kondo 
transparency ^j(w) may be calculated for any pair of 
electrodes (ij). It is a step- wise function in accordance 
with multistage Kondo screening process, but no anoma- 
lous "freezing out" of Kondo effect similar to that shown 
in the left panel of Fig. 21 is expected in this case. 

It is appealing, however, to exploit other specific prop- 
erties of the TQD in the fork configuration. The remark- 
able feature of the mirror asymmetric TQD is that it can 
be viewed as a quantum pendulumi^ii&. This means that 
the superposition of two degenerate states (13) and (23) 
can be considered as a sort of mesoscopic resonating va- 
lence bond (RVB). In these two papers various manifes- 
tations of spin entanglement at even occupation TV = 4, 2 
were investigated. Below we will discuss how this prop- 
erty is manifested in the Kondo regime in the odd occu- 
pation charge sector iV = 1. 

So let us assume that the parameters (gate voltages) 
are tuned so that the TQD is found in a Coulomb block- 
ade valley corresponding to the occupation sector N — 1. 
Within the same approximation as I0J the lowest eigen- 
states I A) of the Hamiltonian Hd ^ for N = 1 are the 
set of spin doublets. 
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(19) 



cos 6* d 



sine '^i- 



V2 



|0), 



with sin 6* — ^/2W/ A ^ 1. The corresponding eigenfunc- 
tions are 



(20) 



Hence, the ground state is a bonding spin doublet Em,, 
and the eigenstates |A) may be interpreted as two even 
and one odd RVB modes of a " spin pendulum" . 

In the low-energy subspace oj <C W'^/A the effec- 
tive spin Hamiltonian which describes the Kondo co- 
tunneling has the same form as (|15|) derived above for 
N = 3. Here, however, all coupling constants are posi- 
tive, 
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and the Kondo temperature is given by Eq. H18|l . 

One may say that in a coherent Kondo-tunneling 
regime the system demonstrates perfect entanglement: 
an electron entering dot 3 from lead 3, splits into two 
components in accordance with the structure of the state 
\Db(j) and this entangled state predetermines the total 
current Ii + I2 through the TQD in the fork geometry. 
Of course, this statement is valid only at zero tempera- 
ture, and one may expect that thermal fluctuations are 
detrimental for a coherent transport, but this effect is 
suppressed as {T/TkY low T. 

The situation becomes even richer, due to the occur- 
rence of soft mode excitations (|20|l : the odd state \Dna) 
may be intermixed with the even state \Dhfj) due to co- 
tunneling process. This intermixing becomes relevant 
provided Tk ^ Eun — Ejja- This inequality is, of course 
invalid for the bare eigenstates (|20|l . but the Haldane 
renormalization of the spectrum similar to that described 
by Eq. (|10|l may result in softening of this mode. It 
might even lead to level crossing provided the tunnel- 
ing rate r„ = PiVi for the non-bonding state \Dna) is 
higher than the rate Tb = piVi cos^ 9 + psV^ sin^ 6 for 
the bonding state \Db<T). If the densities of states are the 
same in all leads, pi = p3, then the condition r„ > Tb 
means Vi > V3. The RG flow trajectories for this case 
are presented in Fig. [S] Like in the cross geometry, the 
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FIG. 5: Haldane flow diagram for the levels Ea of TQD in 
fork geometry, d = nD/Tn- Energy is measured in meV units. 

TQD acquires an additional degeneracy in the critical 
region d ~ dcr- However, in this case the sources of de- 
generacy are the RVB degrees of freedom. As a result. 
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one encounters the problem of Kondo effect due to an 
interplay between spin S and pseudospin '3", where the 
latter describes the pendulum degrees of freedom. This 
problem was discussed in the context of double quantum 
dotSfii triple quantum dotsi2*iS and molecular trimers 
chemisorbed on metallic surfaces^li. In this case the ac- 
tual symmetry of the TQD is SU{A). To show this, we 
derive below the effective spin Hamiltonian following the 
method offered in Ref. 19. 

It is useful to generalize the notion of localized spin op- 
erator S"* — \(T)Ti{a'\ (employing Pauli matrices fi (z — 
x,y, z)) to S\^, = \aA)Ti{a'A'\, in terms of the eigenvec- 
tors \Dba) ,\Dna) from (|19|l . Similar generalization ap- 
plies for the spin operators of the lead electrons: s\y = 

X^kk' "^A kcr''^i'^^''k'o-' ■ Here the index A denotes conduc- 
tion electrons in lead 3 (A = 3), and in leads 1,2 with 
A = e, o corresponding to even and odd combinations of 
conduction electron states in these leads, 



The system of scaling equations has the form: 
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where jb = jbbee, jr, 

jbn = jbnoe- From these equations we derive the Kondo 
temperature 
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'(e,o)kCT ~ y^^'^^''^'^ ^ 4,kcr)- 



(22) 



The pseudospin operators T describing the RVB mode 
are introduced as follows: 

T+ = Y,\Dba){Dna\, = [T+]\ (23) 
r'- = ]^'^{\Db(j){Dha\-\Dna){Dn(j\) . 

a 

The five vector operators Saa' and T constitute the 15 
generators of the SU{A) group. 

Similarly, one may construct the pseudospin operators 
for the electrons in the leads (1,2): 

^+ = r- = [r+]t, (24) 

ka 



The exchange Hamiltonian for TQD with spin and RVB 
degrees of freedom is 



Hsw — E JkXiip^kX ■ + JpT ■ 

k\^P 



(25) 



with K, X — b,n; II, p — e,o, and the coupling constants 
-/kAmp = JxH^t^p = JkXpp are positive Uke in (ED 



Jbbee — 2Jll, Jbb33 — J'u 
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(28) 

with j+ = jb + jn + J3, = jb ~ jn ~_j3- Like in 
the cross geometry, one may manipulate D by chang- 
ing the gate voltages and scan the dependence Tk{D) 
similarly to Fig. |3 This curve has a maximum in 
the critical point and the orbital degrees of free- 

dom are frozen out in the asymptotic regimes D 3> Dcr 
and D <C Dcr- We deal in this case with a symmetry 
crossover SU{2) -> S'C/(4) SU{2). However, there 
is an important difference between the two asymptotic 
SU{2) symmetries. 




FIG. 6: Evolution of Kondo temperature, as determined 
by Eg.ra. 

In the limit d ^ dcr the ground state is Edb, all three 
dots are partially occupied in accordance with the struc- 
ture of the corresponding wave function \Dba) f ll9|l . The 
terms ^ Jbbee , Jbb33 , Jbb3e survivc in the SW Hamiltonian 
(gU, and Eq.(|2Hl for Tk reduces to Eq.ljlHl). In the Umit 
d <^ dcr, the ground state of the TQD is the non-bonding 
state Ejjn with an empty site 3 in accordance with the 
form of the wave function \Dna) H19|) . The SW Hamil- 
tonian (|25|l contains in this case only the term ~ Jnnoo 
and Tk = -Dexp(— 1/jVi)- The Kondo temperature as a 
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function of D has a maximum in a crossing point Den 
but unlike the case of cross geometry (Fig. 13), Tk{D) is 
nonzero on both sides of this maximum (Fig. . 

The change of the ground state wave function from 
the bonding combination |i?bcr) to the non-bonding one 
\Dnij), influences the behavior of tunnel conductance. 
Let us compare the tunnel current between the leads '2' 
and '1' and between the leads '3' and '1', which is defined 
by the components G22 and 6*33 of the three-terminal 
conductance matrix Gij — dli/dVj. 

When calculating these components as a function of 
D, one immediately sees that the Kondo-type ZBA in 
G22 exists in all three regimes, and the peak of this con- 
ductance follows the behavior of Tjf (Z)). The behavior of 
G33 is more peculiar. The ZBA in this channel exists only 
until the even components \Dba) are involved in Kondo 
tunneling. In the limit Tk ^ Eob — Eon at d ^ dcr this 
anomaly disappears, so we encounter a unique situation 
where the Kondo resonance is absent in conductance in 
spite of the presence of Kondo screening. 

However, the resonance Kondo regime in G33 arises as 
a finite bias anomaly (FBA). To describe this tunneling 
one has to retain the terms ~ Jbnoe and Jp in the Hamil- 
tonian (|25|l . These terms describe inelastic tunneling, 
which acquires a form of Kondo resonance at finite bias 
in accordance with the mechanism offered in Ref. I21I 
In that case the tunneling through the double quantum 
dot with even occupation in a singlet ground state was 
considered, and the Kondo regime becomes relevant at 
finite bias, when the difference in the chemical potentials 
of source and drain leads compensates the exchange gap 
between the ground state singlet and excited spin triplet 
state. In our case the FBA arises at odd occupation when 
the bias compensates the energy gap 5 = Ejji, — Ejjn- 
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FIG. 7: Left panel: Tunnel conductance in the channel '2-1' 
for D — Dcr (a), Z) > Dcr (b) and D < Dcr (c). Right panel: 
Tunnel conductance in the channel '3-1' for D > Dcr (a) and 
D < Dcr (b). 

In order to calculate the tunnel conductance in the 
weak coupling regime, we use the modified Golden 
Rule formulas^, which reads for the channel '3-1' as 
G33(ej3,r)/Go ~ |Ji3(ey3,T)|2, where Go = eynh, 
and Ji3 is the solution of the RG fiow equations H27|) 
for D = max{eV3,T}. In case of Kondo-type resonance 
at finite bias eV^, the equation for conductance in the 



vicinity of FBA readaSi 

G/Go - ln"^(max{ey3 - d,T}/TK). (29) 

Similar equation (with V2 — > 0) describes the ZBA in 
tunnel conductance in the channel '2-1'. 

The results of calculation of G22 and G33 are presented 
in Fig. [3 One may see from this figure how the ZBA 
in G22 goes through the maximum at the crossing point 
Dcr whereas the ZBA in G33 transforms into FBA for 
small enough D when S exceeds Tk. 

Discussion of damping effects, which, to a certain de- 
gree, tend to smear the FBA^-'^'^^ is beyond the scope of 
this work. It should be mentioned, however, that it was 
shown in Ref. '2T that there exists a wide enough window 
of parameters, where this damping is not fatal for exis- 
tence of well shaped FBA. It is worth noting also that the 
tunnel conductance between the leads '1' and '2' exists 
in spite of the absence of direct tunneling channel. The 
tunneling mechanism is connected in this case with the 
" pendulum" structure of the electron wave function H19|) 
in TQD. In case of ground state Em this is the bonding 
combination {d\^ + dl^)\0), in case of ground state Eon 
this is the resonating valence bond {d\^ — dl^)\0). In 
the latter case the dot '3' is excluded from cotunneling. 
It is therefore involved only in the determination of the 
Kondo temperature. 



IV. CONCLUDING REMARKS 

We have shown in this paper that triple quantum dots 
in some special geometries demonstrate unusual behav- 
ior in the Kondo tunneling regime. This behavior stems 
from inequivalence of constituents (side valleys and cen- 
tral valley). Asymmetric double quantum dot considered 
in Ref. ^ was the first example of complex quantum dot 
with inner and outer " shell" . Trimers with inequivalent 
side and central valleys studied in Refs. .8, 9, 10. and in the 
present paper are more complicated examples of artificial 
molecules with shell structure. In this case the two side 
dots form an "inner shell" whereas the big central dot 
plays the role of an " outer shell" . If the outer shell is 
open2ii£, it contributes to the indirect exchange between 
the electrons in the inner shell. 

We considered here the case of closed outer shell and 
found that such trimer with odd electron occupation 

= 1,3 possesses properties, which were observed ear- 
lier in dots with even occupation N = 2. In particular, 
the Kondo tunneling may be absent in the ground spin 
doublet state of TQD due to special symmetry proper- 
ties of the wave function (odd I — r symmetry in case 
of TV = 3 and empty outer shell in case of A^ ~ 1). In- 
volvement of the low-lying Kondo-active spin doublet re- 
sults in a two-stage Kondo screening reminiscent of that 
found in quantum dot with occupation N = 2 where the 
spin excitation spectrum is formed by the singlet-triplet 
pair.^'^ Other interesting possibilities now open due to 
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the resonance valence bond structure of the electron wave 
function (|19|l in case of partially occupied inner shell in 
a fork geometry with = 1. In particular the "pendu- 
lum effecfii^ti^ perceived in TQD with even occupation 
TV = 2, 4 may be exploited in this type of TQD as well. 

The transformation of ZBA into FBA under changing 
gate voltage is a special manifestation of general phe- 
nomenon, known as "critical phase transition" , where the 
symmetry of the ground state changes as a function of 



a control parameter. Similar effect should be observed 
in planar and double quantum dots^^ with even occupa- 
tion where the singlet-triplet crossover may occur with 
changing gate voltage or in transition metal molecular 
complexes^^. In the latter case local phonons are essen- 
tially involved in this transitioniSi 
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